Quasi-static and dynamic order-disorder transition in presence of strong pinning 
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We present the results of an experimental study of vortex dynamics in non- twinned Y Ba2Cu^Oa,s,7 
crystal. It is found that critical currents J c and J c .dyn, which correspond to the pinning force in 
the thermal creep and flux flow mode, respectively, non-monotonically vary with the magnetic field. 
However, the minimum in the J c ,d y n(H) dependence is observed in higher fields, compared with the 
minimum position Hod in the J C (H) dependence. Considering that the field Hod corresponds to 
the static order-disorder transition, this difference is explained by partial dynamic ordering of the 
vortex solid. It is concluded that finite transverse barriers guarantee finite density of transverse 
displacements of vortex lines u t ~ c^ao suitable for preservation of the disordered state of the 
moving vortex solid. 

PACS numbers: 74.25.Qt, 74.25.Sv, 74.72.Bk 



The interaction of static and dynamic elastic media 
with chaotic pinning potential is one the chapters of solid 
state physics, which includes dislocations in solids, charge 
density waves, Vigner crystals, and vortex lattices (VL's) 
in Type-II superconductors. The VL's are the most ap- 
propriate objects for the experimental study of elastic 
media, because it is easy to change the strength of pin- 
ning potential in superconductors, as well as the elasticity 
and motion velocity of VL's. An important feature of the 
VL's is the non-monotonous field variation of the pinning 
force F p , which is observed in low-T c (NbSe2 Q, HJ, VjSi 
|) middle-T c (MgB 2 @, Q), and high-T c (BiSrCaCuO 
|, YBaCuO @,[|) superconductors. The increase of the 
pinning force can be explained by softening of the elas- 
tic moduli of VL's in vicinity of the upper critical field 
H C 2 {T) [2| or the melting line H m (T) Q that causes bet- 
ter adaptation of the vortex lines to the pinning land- 
scape. Some alternative models [l(| El suggest forma- 
tion of an ordered vortex solid (VS) in low fields, which 
transforms into a disordered one in some magnetic field 
Hod, though the nature of the order-disorder (OD) tran- 
sition and the mechanism of increasing the force F p may 
be different. These models are supported by correlation 
between the field Hod corresponded to the structural 
OD transition [l| and the onset of the F p increase @ in 
BiCaSrCuO crystals. An actual problem of the VS phase 
is the nature of its ordering under an increased vortex ve- 
locity v. The "shaking temperature" model [l3| suggests 
that transverse vortex displacements u t induced by the 
disorder reduce with increased velocity, Ut oc l/t>; and 
the increase of the velocity above some critical value v c 
results in a dynamic transition from the disordered to or- 
dered state. It was later justified that the increase in 



'Electronic address: Alcksandr.V.Bondarenko@univcr.kharkov.ua 



v leads to a suppression of the pinning in the longitudinal 
(with respect to v) direction only, while pinning barriers 
remain finite in the transverse direction. The effect of 
motion on the transverse barriers, phase state and pin- 
ning force of vortex solid is still controversial issue, and 
this subject first of all requires additional reliable exper- 
imental studies. The goal of this work is experimental 
study of vortex dynamics in the presence of strong pin- 
ning. 

The measurements were performed on detwinned 
YBa2Cu307_,j crystal, annealed in an oxygen atmo- 
sphere at 500° C for one week. Such anneal corresponded 
to an oxygen deficiency 5 ~ 0.13 [l5j and T c ~ 91.8 K. 
The crystal then was held at room temperature for 7 days 
to form clusters of oxygen vacancies, which reduced the 
tension of the field Ho d [IB] • The field variation of the 
pinning force was studied through measurement of the 
current-voltage characteristics, E(J), using the standard 
four-probe method with dc current. The investigated 
sample had rectangular shape with smooth surfaces; its 
dimensions were 3.5x0.4x0.02 mm with the smallest di- 
mension along the c axis; the current was applied along 
the largest dimension; and the distance between the cur- 
rent and potential contacts, and between the potential 
contacts was about 0.5 mm. The measurements were 
performed at a temperature of 86.7 K in the field H || c. 
Fig. [Tj shows the v(J) = cE{J)/B dependencies and 
current variation of the normalized dynamic resistance 
p d (J) = [dE(j)/dJ]/p BS , where p BS = PnB/B c2 $. 
At low currents, the electric field increases exponentially 
with an increase in current and the resistance pd is much 
lower than one. This increase in v and low dynamic resis- 
tance indicate the presence of thermally activated vortex 
creep. At high currents, the v(J) dependence is linear 
and the value of pd is close to 1, indicating the presence 
of the flux flow mode. The critical current in the ther- 
mal creep mode Je can be characterized by the voltage 
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FIG. 1: E( J) curves presented in the linear (a) and semi- 
logarithmical scale (b), and Pd(J) curves presented in the 
semi-logarithmical scale (c) . The inset in panel (b) shows the 
E(J) dependencies measured upon increase (light symbols) 
and decrease (dark symbols) of the current. 



criteria of E = 1 fiV/ cm and E = 100 ^(V/cm, and the 
dynamic critical current J c ,dyn can be determined by ex- 
trapolating the linear parts of the v( J) dependence, cor- 
responded to the flux flow mode, to zero voltage p7| . 
Field variation of the currents Je and J c ,dyn normal- 
ized by their values in a field of 0.5 kOe are shown in 
Fig. [5k.. It is seen that the currents J c and J c ,dyn start 
to increase in the fields above 1.25 kOe and 2.5 kOe, 
respectively, which are substantially smaller in compari- 
son with the fields H C 2 and H m . Therefore this increase 
can not be caused by better adaptation of the vortices to 
the pinning landscape induced by softening of the elas- 
tic moduli. Obtained field variation of the currents Je 
and J c ^dym and the peculiarities of vortex dynamics can 
be explained in frames of the model proposed by Ertas 
and Nelson [ioj |. It is assumed that the OD transition 
occurs when transverse displacements of vortex lines ex- 
ceed the value of c L ao, where ao — \/$o/-B is inter- 
vortex distance, $o is the flux quantum, and cl is the 
Lindemann number. The field is defined by equality of 
energies E c i(Hod) — E p {Hqd), where E p is the pinning 
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FIG. 2: (a) Field variation of the current J c ,dyn and Je 
normalized by their values in a field of 0.5 kOe. (b) Field 
variation of the velocities v p and v min correspondent to the 
peak and minimum position in the pd{J) dependencies, re- 
spectively. The inset in panel (b) shows sketch of the trans- 
verse vortex displacement u t ,L correspondent to the Linde- 
mann criteria. Dash and solid circles correspond to the lower 
and upper boundaries of the displacements ut,L (see the text), 
respectively, in the static VS in magnetic field H > Hod- 
Dot ellipses show evolution of the maximal displacements Ut,L 
upon increase of the velocity v. Dashed region corresponds 
to the displacements u t ,L in the dynamic VS. 



energy, E e i ~ c L seoao is increase of the elastic energy 
caused by displacements Ut — clclq, e is the anisotropy 
parameter, Eq = (^q/AttX) 2 is the line tension of vortex 
line and A is the penetration depth. As evident from 
Fig. [3Ja, the minimum position does not depend on the 
driving force within the creep regime in agreement with 
magnetization measurements 0, H[ ■ This means that the 
value of ratio E e i / E pi and, therefore the energy E p , is not 
changed, indicating that minimum in the Je(H) curve 
corresponds to static OD transition, Hod — 1-25 kOe. 
Estimations presented below show that static OD tran- 
sition in our sample is caused by vortex interaction with 
the clusters of oxygen vacancies rather than with the iso- 
lated oxygen vacancies. Indeed, for the point disorder the 
pinning energy is 0, El E p ~ ( 7 £ 2 £ £ 4 ) 1/3 (W£c) 1/5 , 
where L ~ 2ea is the length of longitudinal fluc- 
tuations, L c ~ ££,{Jo/ Jd) 1 / 2 is the correlation length, 
Jo = 4ce /3\/3£$o ^ s depairing current, £ is the co- 
herence length, and 7 ~ (J c &o/c) 2 L c is the disorder pa- 
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rameter. Using realistic for the YBa2Cu307_5 supercon- 
ductor parameters (A = 500 nm, £ = 4 nm, and e = 1/7) 
and experimental value of the depinning current J c ,d y n < 
5 kA/cm 2 we obtain the energy E p < 2-10 -16 erg, which 
is about 25 times smaller compared to the elastic energy 
E e i — c^eeoao ~ 5-10~ 15 erg estimated for the cl = 0.2 
and Hod — 1-25 kOe. The pinning energy induced by 
vortex interaction with the clusters of oxygen vacancies 
equals the condensation energy U c w (H^ / '8tt)V c i , where 
H c = <f>o/2v / 27rA£ is the thermodynamic critical field and 
V c i is the volume of clusters. For spherical clusters with 
radius r~(we obtain the energy E p ~ U c « 10~ 14 erg, 
which is suitable for occurrence of the OD transition in 
the field of 1.25 kOe. 

As it is shown in Fig. [2b, minimum in the J c ,dyn(H) 
curve occurs in a field of 2.5 kOe, which is about two 
times exceeds the value of Hod- Also, above the mini- 
mum position, the current J c ,dyn(H) increases with the 
field more gradually in comparison with increase of the 
current Je above the OD transition. This difference can 
be explained by suppression of the longitudinal and con- 
servation of the transverse pinning barriers, as it was 
theoretically predicted in In frames of the "shak- 

ing temperature" model [13j , this means conservation of 
the transverse (with respect to vector v) u± and reduc- 
tion of the parallel u\\ oc 1/v component of the displace- 
ments Ut = (w||) 2 + (uj_) 2 with increased velocity v. In 
magnetic field H || c and in presence of the chaotic pin- 
ning potential, spatial distribution of the displacements 
is isotropic; and in the field H > Hod, the displace- 
ments Ut.L, which correspond to the Lindemann criteria, 
fall in the interval clcio(Hod) > Ut.L > ClCLo(H), as it is 
shown schematically in the inset of Fig. 2b. Density of 
the displacements (the number of vortex displacements 
Ut,L P er unit length of vortex line) n t% L is proportional to 
the area of ring confined by the upper (solid circle) and 
lower (dashed circle) boundary of the displacements u^l- 
Reduction of the component u\\ with increased velocity 
v leads to reduction of the upper boundary (dotted lines 
for velocities v-i > v\ ^ 0) , and thus to reduction of the 
density n t ,L- It is important, that for any finite velocity 
v the component uu is finite, and thus the cross-hatched 
area at the diagram, which corresponds to the displace- 
ments u^l, and the density n t £ is also finite. Increase of 
the field reduces the lower boundary of the displacements 
Ut.L, and therefore the density n t .L increases. 

Considering that the displacements Ut.L produce the 
dislocations in the VS phase, and increase of the density 
n t .L results in an increase of the current J c ,dyn 0, the 
field variation of the currents Je and J c ,dyn can be ex- 
plained in the following way. In low fields, the ordered 
VS phase, which is characterized by the absence of dis- 
location and realization of the ID pinning, is formed, 
and the currents Je and J c ,d y n decrease with increased 
field due to enhancement of the vortex- vortex interaction, 
making difficult to fit the vortices in the pinning land- 
scape. Above the OD transition, the VS phase contains 
dislocations that results in occurance of the 3D pinning 



[l(|, and thus the current Je increases at the transition 
point Hod due to dimensional crossover in the pinning 
[20l l2l| . Further increase of the current Je with mag- 
netic field is caused by increase of the density n^L, as it 
was found in (l9| . The density nt : L is smaller in the mov- 
ing VS phase than in the static VS phase, but it is finite 
and increases with the field. Therefore, the J c .dyn(H) 
dependence is determined by competition between de- 
crease of the pinning force caused by enhancement of 
the vortex-vortex interaction and increase of the pinning 
force associated with increase of the density nt^L- In our 
measurements, the former mechanism dominates in mag- 
netic fields H < 2 kOe, while the last one dominates in 
the fields H > 3 kOe. 

Proposed in terp retation agrees with numerical simu- 
lations of 2D HI 11 H Hi and 3D [H VL's in the 
presence of strong pinning. First, it was shown that in 
the flux flow mode the disordered state of the VL's is 
preserved [l5 L 22^ 24|, [25| , and the transverse barriers re- 
main finite |23l. |24(. Second, the v(J) curves cross one 
another near the OD transition [15| . Third, our inter- 
pretation implies that cross-hatched area in the diagram 
collapses to a segment at v — > oo, indicating that mov- 
ing VS can be ordered in agreement with conclusion in 
[Fa ]. Finally, the onset of ordering of the moving VS 
phase is manifested as a peak in the Pd{J) curves, and 
the end of ordering corresponds to value of the resistance 
Pd{J)=l [12, Hll, and in our measurements peak in the 
Pd(J) curves appears in the fields H > Hod- Following 
computer simulations, we determined the field variation 
of the velocities v p and v m in, which correspond to the 
peak and minimum positions in the Pd{J) curves respec- 
tively. As it is shown in Fig. 3b, the velocity v p and the 
difference Av = v m m — v p increase with the field indicat- 
ing that the critical velocity of the ordering as well as the 
interval of velocities Av, in which the ordering realizes, 
increase with the field. This behavior is plausible consid- 
ering that the lower boundary of the displacements Ut t L 
decreases with the increased field that requires higher v's 
to decrease the amplitude below this boundary. Also, the 
difference between the upper and lower boundary of the 
displacements Ut,L, Au = CL[ao(-ffozj)~ a o(^)], increases 
with the field that results in increase of the difference Av. 

Our interpretation allows explaining occurrence of the 
hysteresis effect in the curve v{ J) measured with the in- 
creased and decreased current in a field of 1.5 kOe, and 
absence of the hysteresis effect below and quite above 
the OD transition. Indeed, in close vicinity to the OD 
transition, (HjHoD — 1) << 1, the density n t .L in the 
dynamic VS is much smaller than in static VS, and small 
increase of the velocity v leads to dynamic transition into 
the ordered state. In this case the "shaking tempera- 
ture" model predicts the hysteresis effect, which reflects 
the "overheated state" of the ordered dynamic VS. The 
decrease in density ti^l quite above the OD transition 
is not dramatic, and transition from strongly disordered 
static VS to less disordered dynamic VS occurs in a wide 
interval of velocities Aw without hysteresis. It is impor- 
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tant to notice that the E(J) curves measured after zero 
field cooling coincide with the E( J) curves measured af- 
ter non zero field cooling, that indicates the absence of 
metastable states in the VS. This agrees with experi- 
mental studies of the YBaCuO crystals: the metastable 
states exist in vicinity of the vortex sold - vortex liquid 
transition, but they disappear below this transition [261 ] . 

Recent quantitative theory of the dynamic VS by 
Rosenstein and Zhuravlev [ll[ predicts jump- like increase 
of the pinning force at the OD transition. It is evident 
that this theory does not describe our results because 
increase of the currents Je and J c ,dyn(H) occurs in dif- 
ferent fields, and field variation of the currents does not 
show the jump-like increase. 

The obtained field variation of the currents Je and 
Jc,dyn, occurrence of the hysteresis effect in close vicin- 
ity to the OD transition, and absence of the metastable 
states in the VS are different from that in superconduc- 
tors with weak bulk pinning. For example, in crystals 
NbSe2 [23 and MgB 2 Q, the current J r ..dyn increases 
in a jump- like manner at the OD transition [23, HH, an d 
the hysteresis effect occurs in rather wide interval of mag- 
netic fields and it is caused by presence of the metastable 



states in the VS |5L|27|], which are induced the effect of 
surface barriers 2{|. The surface barriers in the NbSe2 
[30l | and MgB2 [I | cause asymmetry of the magnetization 
loops, and this asymmetry reflects a difference in the bar- 
riers for vortex entrance and exit of samples [3lJ. The 
magnetization loops of the YBaCuO crystals are sym- 
metric indicating negligible effect of the surface barriers. 
Therefore obtained field variation of the current Je cor- 
responds to equilibrium quasistatic VS. 

In conclusion, we determined field variation of the crit- 
ical currents in the quasistatic and dynamic vortex solid. 
The currents non-monotonously vary with the field, but 
minimum position in the J c .d y n(H) dependence is shifted 
to higher fields in comparison with the minimum in the 
Je{H) dependence. The difference is interpreted by par- 
tial ordering of the vortex solid with increased vortex ve- 
locity. The disordered state of the dynamic vortex solid 
is attributed to preservation of finite transverse pinning 
barriers that guarantees presence of the transverse vor- 
tex displacements suitable for formation of dislocations. 
This interpretation allows explaining observed increase 
of the critical current of the dynamic ordering. 
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